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ABSTRACT 

In recent years, the concept of using the rudder for roll stabilization of a ship in 
a seaway has been investigated with good results. Such designs have been used to 
solve the roll problem of a ship on steady course while compensating for yaw 
perturbations. 

To complete a review of the possible design strategies to meet the specifications 
for this model, the effects of the feedback gains on rolling and yawing are studied in 
detail. Roll angle and roll rate feedback are used to control the rudder. 

Roll stabilization with the rudder in various sinusoidal sea states is studied by 
simulation on the IBM digital computer. The model used is based on the data 
obtained from a typical naval ship. The Root Locus method is used to design the 
feedback gains. The computer simulation programs are written in Digital Simulation 
Language (DSL/VS), are plotted as data in DISSPLA and include the effects of rudder 
servo nonlinearities, which seriously restrict the ability of the rudder to reduce roll. 
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I. INTRODUCTION 



There is a need to reduce roll on various types of ships in order to increase the 
comfort of the passengers, to provide safe operation of container ships in bad weather 
and to provide stable platforms for helicopter operations. 

Passive or active tanks, or stabilization fins are most commonly used at present 
to control roll. But both have disadvantages: tanks require a lot of space, fins 
introduce a considerable drag, are expensive, and also require space for the hydraulic 
actuators. 

Recently there has been a considerable amount of interest in Roll Stabilization 
systems which use the rudder (RRS), since rudder motions not only affect a ship's 
heading but influence the roll motions as well. This resulted in the successful 
installation of a rudder roll stabilization system on some ships of the United States 
Coast Guard Cutter HAMILTON class (378 foot) as reported by references (1) and 
( 2 ). 

This thesis describes a simple mathematical model for the transfer between the 
rudder angle and the two outputs : Rate of turn (heading rate) and roll angle. 
Computer simulations demonstrate results of a design procedure for a combined 
controller for roll stabilization with the rudder. Models are also provided for the 
disturbances and the steering machine. 

This thesis is organized as follows : Chapter II gives mathematical models of the 
ship, the disturbances and the steering machine. Chapter III deals with the controller 
design. Chapter IV gives results of computer simulation and Chapter V summarizes the 
results giving conclusions and recommendations for further research. 
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II. MATHEMATICAL MODELING 



The motions of a ship in waves depend on 

• The dynamics of the ship 

• The disturbances and 

• The controller output, which is influenced by the steering machine. 

A. THE SHIP'S DYNAMICS 

In van Amerongen and Van Cappelle, 1981 [Ref. 3: p. F2 1-3], the basic 
mathematical model was derived which describes the dynamics between the rudder as 
input signal and yaw and roll as output signals. The results will be summarized below. 

The basic equation are : 



Y = m (v-ur). 



( 2 . 1 ) 



K= I x O. 



( 2 . 2 ) 



N = I z r. 



(2.3) 



where 



m 



is the mass of the ship, including the mass of the 
displaced water. 



U=u. + v. + w. : Linear velocity vector with components along x, y, 

1 ] K 



and z-axis. 



w 



u 



V 



is rate of surging, 
is rate of swaying, 
is rate of heaving. 
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12 = p. + q. + r k : Vector angular velocity with components about 
x, y, and z-axis respectively. 



p is rate of roll ( = <P). 

q is rate of pitch ( = O). 

r is rate of yaw ( = T). 

F=x. 4-v. 4-z. : Vector force acting on the ship. 

1 " J * 

x is hydrodynamic force along x-axis. 

v is hydrodynamic force along y-axis. 

z is hydrodynamic force along z-axis. 

M= K. + M. + N. : Vector moment acting on the ship. 
1 J K 



K is rolling moment about x-axis. 

M is pitching moment about y-axis 

N is yawing moment about z-axis. 

Angular moment = I x pj + 1 q. + I z r fc . 



I is mass moment of inertia about x-axis. 

A 

I y is mass moment of inertia about y-axis. 

I is mass moment of inertia about z-axis. 

Z 



The coordinate system where the above variables are defined is shown in Figure 2.1. 
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Figure 2.1 Coordinate System of the Ship. 
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The left-hand sides of equation (2.1) - (2.3) can be expanded into a Taylor series 
in which the higher - order terms are neglected : 

Y = Y y v + Y f r + Y^O + Y 6 6 (2.4) 

K = K v v + K f r + K^O) + + K^S (2.5) 



N = N y v + N f r + (2.6) 

Laplace transform of equations (2.1) - (2.6), substitution of equations (2.1) and 
(2.4) in (2.5) and (2.6) and substitution of equations (2.2) and (2.3) in (2.5) and (2.6) 
yield, under the assumption that v is small : 

<l> < s > = ^,,”"1 2 (K s 6(5) - K f r(s)}. (2.7) 

s i + 2C(0 s + o> i v ' 

^ n n 

r(s) = {N 6 8 (s) - O(s)}. (2.8) 

ST f +l ° v 

These equations can be illustrated with the block diagrams of Figure 2.2 and 
Figure 2.3. 

The model used in the subsequent analysis was "derived" as follows. First it was 
assumed that the influence of the rudder on heading rate would not be strongly 
influenced by roll angle. The conclusion of this assumption is that a heading rate 
model could be established using only rudder angle as input. The second assumption 
was that linear models would be sufficient to identify the major features of the 
dynamics that would be required. With the further assumptions that the roll dynamics 
could be modeled as a second order system with inputs of rudder and heading rate 
only, the trial model looks like equation (2.7) and (2.8). 
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Figure 2.2 Block Diagram for Ship and Control System 
(With Compensator). 



where 

<P is rolling angle, 

r. is heading rate. 

6 is rudder angle. 

N s is a constant relating rudder angle to heading rate. 

is a constant relating rudder angle to rolling. 

K r is a constant relating heading rate to rolling. 

is a constant relating rolling rate to heading rate. 

T f is the time constant of the heading rate system. 

S. is the Laplace transform variable ( = d/dt). 

C, is the damping ratio of the roll dynamics. 

co n is the natural frequency of the roll dynamics. 

The parameter values of this model have been estimated from full-scale trials 
with a naval ship. Table 1 gives some results [Ref. 4: p. 44 1. Note that N(j) has been 
assumed to be zero. 
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Figure 2.3 Block Diagram of the Dynamics between Rudder and 
Yaw andTloll (Without Compensator). 

TABLE 1 

PARAMETERS FOR A NAVAL SHIP 





at 13 knots 


at 17.5 knots 


at 20.5 knots 


n 5 


0.050 


0.077 


0.0S7 


k 6 


0.13 


0.20 


0.25 


K 

r 


4.2 


4.9 


5.3 


T r 


6.8 


6.7 


5.9 


5 


0.18 


0.18 


0.22 


( °n 


0.59 


0.59 


0.58 



(Note : All the above coefficients assume that the basic unit of time is the second. The 
basic unit of angular displacement can be either radians or degrees without altering the 
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coefficients. However, for convenience, in the remainder of this thesis all angles will be 
in degrees.) 

It is interesting to compare the coefficients obtained for the naval ship in Table 1 
[Ref. 4: p. 44], with those of a ship for which a Rudder-Roll-Stabilization has 
apparently been successfully implemented. The comparison data contained in Table 2 
was obtained for the CONFIDENCE from John R. Ware [Ref. 1: p.ll]. There is a 
remarkable similarity in most of the parameters, there is a slight difference in the 
coefficient relating rudder angle to turn rate. Clearly, the CONFIDENCE will turn 
much more rapidly, with a consequent reversal of the rudder induced roll. 

TABLE 2 

PARAMETERS FOR CONFIDENCE MODEL 





at 16 knots 


N 


0.175 


K s 


0.17 


K r 


4.57 


T r 


8.0 


$ 


0.195 


<°n 


0.565 



B. THE DISTURBANCES 

Waves are the most important disturbance with respect to roll. They can be 
described by means of a frequency spectrum, for instance the Bretschneider Spectrum 
[Ref. 5], or the Pierson-Moskowitz Spectrum [Ref. 6]. However, these models of the 
sea state are not conveniently used for the design of a roll stabilizing compensator. In 
this thesis the wave disturbance used is a single frequency sinusoid, with co = 0.5 
rad/sec which is the natural frequency of the ship in roll. This disturbance is used to 
cause roll motion, and the compensator is adjusted to provide maximum reduction in 
roll amplitude. 

Two methods are used to verify the behavior of the system over the range of 
frequencies in a normal sea state. These were 
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(1) A number of single frequency sinusoids, of different frequencies and properly 
adjusted amplitudes were used individually to determine the spectral response 
as in Figures 4.21 through 4.25. 

(2) To approximate a real sea, eight sine waves of appropriate frequencies and 
amplitudes were added and the resulting wave used to check the roll 
suppression, as shown in Figures 4.3 through 4.12. 

C. THE STEERING MACHINE 

In order to reduce the roll motion of a ship using its rudder, the steering system 
must meet certain economic, hydrodynamic and machinery criteria. These criteria can 
be met by using as much of the existing expensive steering system machinery as 
possible, thus voiding most of the capital and maintenance costs associated with 
conventional roll stabilizers. 

When the rudder is going to be used for reduction of roll motions it should be 
able to follow frequencies near the natural roll frequency, o> , without- a noticeable 
phase lag. The steering machine used can be described by the simplified block diagram 
of Figure 2.4a, and autopilot can be described by the block diagram of Figure 2.4b. 




Figure 2.4a Block Diagram of the Steering Machine 
KRSER = 10.0. 

A positioning system is used to position the rudder for small signal operation. 
Its bandwidth should be much larger than that of the ship. In order to prevent a 
phase lag, the angle of rotation and the angle rate have limits. 
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The rudder angle limit is determined either by the mechanical construction of the 
steering machine or is set by the autopilot. The rudder rate limit is determined by the 
construction of the steering machine and by the number of hydraulic pumps which are 
in operation. 

In the design of the controller care must be taken to ensure that the derivative of 
the output signal of the controller is less than the maximum rudder rate in order to 
prevent phase lag. 
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III. CONTROLLER DESIGN 



A. INTRODUCTION 

The primary objective of this chapter is to present procedures for the design and 
compensation of the wave input to roll output control system. Compensation is the 
adjustment of the system in order to limit the maximum roll magnitude. 

The approaches to control system analysis and compensation used in this thesis 
arc the Root-locus approach and the frequency-response (Bode) approach. The system 
was tested in 2 phases. In the first phase, the uncompensated system is tested by 
adjusting the frequency of the wave input. The wave function was generated as a 
sinusoidal input. Since both roll angle and roll rate are easily measured on ship, it was 
proposed that these be fed back to provide roll damping. 

The system which may be designed by a trial-and-error approach is checked to 
see whether the designed system satisfies a desired maximum roll magnitude. 

B. ROOT-LOCUS APPROACH TO CONTROL SYSTEM DESIGN 

The Root- Locus method is a graphical method for determining the locations of 
all closed loop poles from knowledge of the location of the open loop poles and zeros 
as some parameter (usually the gain) is varied from zero to infinity. 




Figure 3.1 Block Diagram Manipulation for Root Locus Method 
(Wave Input to Roll Output). 
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In Figure 3.1 



0) = 0.59 

n 

C = 0.18 

KRSER = 10.0 

In Figure 3.1, the transfer function of the roll compensator, G c , should be 

G c (s)= K 2 *S+ Kj = K 2 (s+ Kj/K 2 ). (3.1) 

In order to solve this problem, the Root-locus method was used to check the 
feasibility of this proposal. First, the ship's dynamics are expressed as a rolling transfer 
function. That is, the feedback term is neglected. 

In Figure 3.1, the roll to wave transfer function, G $ , is 

co 2 0.3481 

G,(s) = -5 0 =- = — 5 (3.2) 

s s 2 + 2^(0 s + co 2 s 2 -4- 0.2 1 24s T- 0. 348 1 

Also, we can get the Root Locus for the roll transfer function. The result is 
Figure 3.2. 

The rudder servo is very fast compared to the natural frequency of the ship's 
rolling mode, and the transfer function of the rudder servo, G f , may be approximated 
as 



G r (s) = 



KRSER 
S+ KRSER 



10 

S+ 10 



The loop transfer function, Gj , is 



G lo( S ) 



3.481(K 2 s+Kj) 

(s+ 10)(s 2 + 0.2 124s + 0.3481) 



(3.3) 



(3.4) 



Note that the open loop has three poles and one zero. The poles are at fixed 
locations, and the location of the zero can be chosen as desired. This leads to a Root- 
Locus as shown on Figure 3.3, where the location of the zero has been chosen 
arbitrarily to illustrate the behavior of the roots. The actual location of this zero is 
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determined by the coefficients Kj and K.,. Root locations are changed by setting the 
gain K r Proper choice of and K-, determines the £ and to n of the roots and the 
damping of these roots is, of course, the damping of the ship roll motion. 

The closed loop transfer function, G cl , is 

0.3481(s + 10) 

ci S ” s 3 + 10.2 124s 2 + (2.4721 + 3.48 1 K 2 )s + 3.481(1 + Kj) 

From equation 3.5: 

The zero at S=-10 should have little or no effect, since it is outside the 
bandwidth. Therefore the response depends on the roots of the denominator and on 
the D.C gain. 

From equation 3.4: 

When we feedback roll 4- roll rate, we change the location of the zero. So the 
range of the zero is from ”0" to "10.0" in calculation. For a small zero the root locus 
will have a loop near the origin, will go to the real axis, then out to the asymptote at 
about -5.0. When the zero is out near -1.0, the Root-Locus from the complex poles 
will not go to the real axis, but will stay complex as it goes to the asymptote. 

The best we can do is use the rudder to improve the damping of the roll. Then 
small zeros are not useful because we do not want to change co n . 
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ROOT LOCUS PLOT 



WITHOUT COMPENSATOR 




Figure 3.2 .Root Locus for Wave to Roll 
(Without compensator). 
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ROOT LOCUS PLOT 



WITH COMPENSATOR 




Figure 3.3a Root Locus for Wave to Roll 
(With Compensator) 
z = 1.0, Kj = 4 . 0 , K, = 4.0. 
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ROOT LOCUS PLOT 



WITH COMPENSATOR 




Figure 3.3b 



Root Locus for Stabilized Model (Main Part) 
_ (With Compensator) ^ ^ 
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C. FREQUENCY-RESPONSE (BODE) APPROACH TO CONTROL SYSTEM 
DESIGN 

In dealing with the problem of compensating control systems via frequency- 
domain techniques, we control the transient-response behavior in terms of such 
frequency-domain specifications as phase margin, gain margin, resonant peak value, 
and bandwidth. Design in the frequency domain is indirect because the system is 
designed to satisfy these frequency-domain specifications rather than time-domain 
specifications. 




Figure 3.4 Block Diagram Manipulation for Bode Method 
(Wave Fnput to Roll Output). 



In Figure 3.4 

0) = 0.59 

n 

5 = 0.18 
KRSER = 10.0 

From the block diagram on Figure 3.4, the closed loop contains a rudder servo, 
and a compensator block to stabilize the system and provide damping of system 
response to the wave input. The loop transfer function, G, o , is 




K,((K 2 /K,)S+1) 

(0.1S+ 1)(2.8727S 2 + 0.6102S + 1) 



(3.6) 
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If Kj > 1, then Bode gain increases and the gain crossover changes accordingly. 
The transfer function, equation 3.6, for G ]o (s) can be written as 



S( KyK^+l 



G. (s) = [ K, ( — - — - — — }( }] (3.7) 

10 1 (S/10)+l (S 2 / 0.3481)+(0.2124/0.3481)S+ 1 



That is, the feedback (G c ) and rudder servo (G f ) transfer functions can be 
combined to have the algebraic form of a lead filter. Also, we can get the frequency 
responses (Bode diagram). The results are Figure 3.5 to Figure 3.7. Figure 3.5 shows 
the uncompensated roll to wave frequency response and Figure 3.6 and Figure 3.7 
show the roll to wave frequency responses (open loop and closed loop) with 
compensator (Kj=4.0, K 2 = 4.0). From this, it was observed that the ship has a 
natural frequency of in roll 0.59 radians per second which corresponds to a period of 
10.65 seconds. 

Then the effect (seen on the Bode diagram in Figure 3.6) of any gain, Kj, is to 
raise or lower the |G| curve. Raising the |G| moves the gain crossover to higher 
frequency, regardless of the location of the zero of the filter. 

The purpose of the filter zero is to increase damping of the roll frequency. To do 
this we want to introduce positive phase at the gain crossover. If we put the zero at a 
frequency lower than the gain crossover (caused by the gain Kj), we increase the 
bandwidth which is undesirable. 

However, if we place the zero at a higher frequency, we increase the damping 
(phase margin) by a reasonable amount without appreciable increase in bandwidth. 

These results are essentially the same as shown by Root-Locus analysis. 

In Figure 3.6, the pole of the rudder servo becomes a zero of the closed loop, and 
the effect of compensation is to move the real root and complex roots to new 
locations, presumably with larger negative real parts for the complex roots. If the 
increase in gain has been small, on the root locus the roots have not moved very far. 

Therefore the real effect of the compensator has been to move the complex roots. 
The effect can be seen on the closed loop Bode diagram. 

The above implies that roll would be reduced at low frequency, though not very 
much, but would be reduced substantially in the range of frequencies near resonance. 
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Figure 3.6 Frequency Response in Open Loop for Wave to Roll 

(With Compensator) 

Z = 1.0 ,Kj = 4.0, K 2 = 4.0. 



31 




CLOSED LOOP 
BODE PLOT 
WITH COMPENSATOR 



(»»P) ZSVHd 




(ap) aanjJNOYPi 



Figure 3.7 Frequency Response in Closed Loop for Wave to Roll 

' (With Compensator) 

Z = 1.0, K t = 4.0, K 2 = 4.0. 
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FREQUENCY (cj, rmd>ec) 



IV. THE RESULTS OF COMPUTER SIMULATION 



A. SIMULATION TEST STUDIES 

This chapter will give a few results of the computer simulations. During the 
computer simulation, the parameter values used in this thesis are taken from in 
Table 1 with ship speed at 17.5 knots. The autopilot gain constant and rudder servo 
sain constant are given in Fisure 2.4. 

The wave disturbance was generated as a sinusoidal input which produced a roll 
angle of ± 20 degrees in the open system at the resonant frequency which is 0.5 
rad/sec. This was used for analysis and design. A complex wave consisting of a sum 
of sinusoids of different frequencies and amplitudes (see program, Appendixces E and 
F) was used to test the chosen designs. The test waves are shown in Figure 4.1. These 
input waves were used to drive the 'open system', the system with ideal rudder, and the 
system with real rudder. 

Figure 4.2 used the single sine wave to test the system with ideal rudder. 
Response of the compensated system is compared with that of the uncompensated 
system. In the ideal system with compensator in this reduces the maximum roll angle 
to approximately ± 9 degrees. So the roll reduction is 55 percent. The numerical 
results are given in table 5 in Appendix I. System response to the complex wave is 
given in Figure 4.3. 

It appears that a considerable reduction with ideal rudder (approximately 71 
percent : Z = 2.0, Kj = 8.0, K, = 4.0) can be obtained, while the heading deviation is 
increased a small amount. However, it appears that a rudder angle of ± 32.6 degree 
and rudder speeds of ± 16.3 deg/sec are required to achieve this. 

In practice there are limitations on both the maximum rudder angle and the 
maximum rudder speed. The maximum rudder angle limits the maximum moment 
which can be applied. Figures 4.4 through 4.8 compare the behavior of the open 
system with that of the rudder roll stabilized ship when the rudder angle is limited to 
± 20 degrees and the rudder rate limit is ± 10 deg/ sec. Each Figure shows the effect 
of a different compensation design. The system of Figure 4.5 appears to have 
minimum rudder activity and minimum yawing, but (from Table 5, Appendix I) the 
maximum roll angle is slightly larger than for the other design. 
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Data for various real rudder rate limits and various rudder angle limits were 
obtained and are given in Table 5. From this table, we can summarize the results as 
shown in Table 3 : 



TABLE 3 

MINIMUM VALUES OF <J> NIAX WITH DIFFERENT ZERO LOCATIONS 

O) = 0.5 



U.«< de &' sec > 


6 (dee) 

max' 


z 


. K 1 _ 


K, 


T 

max 


<J> 

max 


5.0 


5.0 


1.0 


4.0 


4.0 


0.2 


17.9 


10.0 


1.0 


4.0 


4.0 


0.5 


16.2 


20.0 


1.0 


4.0 


4.0 


0.6 


17.4 


30.0 


1.0 


4.0 


4.0 


0.8 


17.4 


10.0 


5.0 


1.0 


4.0 


4.0 


0.3 


17.7 


10.0 


1.0 


4.0 


4.0 


0.6 


14.9 


20.0 


1.0 


4.0 


4.0 


1.0 


10.5 


30.0 


1.0 


4.0 


4.0 


1.2 


9.0 


15.0 


5.0 


1.0 


4.0 


4.0 


0.3 


17.7 


10.0 


1.0 


4.0 


4.0 


0.6 


14.8 


20.0 


1.0 


4.0 


4.0 


1.1 


10.5 


30.0 


1.0 


4.0 


4.0 


1.4 


9.0 


20.0 


5.0 


1.0 


4.0 


4.0 


0.3 


17.6 


10.0 


1.0 


4.0 


4.0 


0.6 


15.0 


20.0 


1.0 


4.0 


4.0 


0.6 


10.5 


30.0 


1.0 


4.0 


4.0 


0.6 


9.0 



It is clear that considerable reduction in roll angle can be achieved. A major 
requirement is high rudder rate. However for a given rudder rate good roll reduction 
can be obtained with several locations of the compensator zero. It may also be 
observed that roll reduction is accompanied by increased yawing unless the rudder rate 
is quite high. Thus a number of options are available. 

Maximum rudder rate is determined by the design of the rudder actuator system, 
but rudder rate less than the maximum could be set in the autopilot. In like manner 
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rudder angle limits can be set in the autopilot, and with proper design the feedback 
gains Kj and K-> can be made adjustable. Thus it is possible to design a rudder roll 
stabilization system which has considerable capability for roll reduction, but permits 
the user to choose the specific operating condition. 

Figures 4.9 through 4.12, shows system response for different feedback gains 
and K 2 ) with the same zero location when a complex wave causes rolling. Rudder 
activity and yaw motion are compared with that which would result if no stabilization 
was used. 

Figures 4.13 through 4.16 show the effect of rudder rate on the maximum roll 
angle when the maximum rudder angle is limited. It is clear that for small rudder angle 
limits there is very little reduction in roll angle, and the rudder rate and zero locations 
have no noticeable effect. As the rudder angle limit is increased greater roll reduction 
becomes available providing the rudder rate is increased. From Figure 4.16 it is seen 
that maximum roll reduction is obtained when both rudder angle and rudder rate are 
maximum. Under these conditions the location of the compensator zero is important 
and from the data available it appears that the zero locations should be in the range of 
1.5 < Z < 2.0. 

Figures 4.17 through 4.20 give a simulation study of the effects of roll 
stabilization on ship yawing. It is seen that those conditions which provide maximum 
roll reduction also result in increased yawing. Thus the user must choose operating 
values which best suit his needs in a particular sea state. 

Design of the roll stabilization was accomplished by considering rolling motions 
at the natural frequency only. To verify the results, the system was tested at a number 
of frequencies within the wave spectrum. Results are illustrated in Figures 4.21 
through 4.25. Rolling is maximum when the system is 'open system' ie, no 
stabilization. Suppression of roll is greatest with an ideal rudder. When the rudder 
angle and rate are limited the effectiveness of the stabilization is reduced. A high 
rudder rate is essential, and if this is combined with a large rudder angle, roll reduction 
is significant over the entire spectrum. 
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B. SIMULATION RESULTS 




Fieure 4.1 Performance of Roll Angle for Open System, 
Sine and Complex Wave Input. 
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Performance of Ideal System, Compensated, 
Sine Wave Input 

K.j-4.0, K.,-4.0, (o = 0.5 rad/sec. 
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OPEN SYSTEM BNO IDEAL SYSTEM' SINE WAVE INPUT 
K1-4.0. K2-4.0. W-0.5 RAD/SEC 
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Figure 4.3 Performance of Ideal System, Compensated, 
Complex Wave Input 
Z= l.S, K,=4.0, K 2 = 4.0. 
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Figure 4.4 Performance of Real System, Compensated, Complex Wave Input 

Z = I.S, k, = 4.5, k., = 2.5 

Maximum Rudder Rate: ± 10 deg, sec, Maximum Rudder Angle : ±20 deg. 
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Figure 4.5 
Maximum Rudder Rate 



Performance of Real System, Compensated, Complex Wave Input 
' = 0.5, K 1 = 2.0.K, = 4.0 



Z v,.^, ,v 1V -,. v 

± 10 degree, Maximum Rudder Angle : ± 20 deg. 
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Fiaure 4.6 Performance of Real Svstem, Compensated, Complex Wave Input 

Z= 1.0, K. =4.0. K, = 4.0 

Maximum Rudder Rate : ± 10 deg, sec, Maximum Rudder Angle : ±20 deg. 
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Figure 4.7 Performance of Real System, Compensated, Complex Wave Input 
Nf . Z= 1.5, K, = 6.0. K, = 4.0 1 F 

Maximum Rudder Rate : ± 10 deg, 'Sec, Maximum Rudder Angle : ±20 deg. 
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Figure 4.8 Performance of Real System, Compensated, Complex Wave Input 

Z= 2.0, K. = 8.0 K, = 4.0 

Maximum Rudder Rate : ± 10 dcg/$ec, Maximum Rudder Angle : ±20 deg. 
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Figure 4.9 Performance of Real System, Compensated, Complex Wave Input 

*** 1 J Fv i 1 . 0 , k ^ 1 .0 

Maximum Rudder Rate : ± 10 deg/ Sec, Maximum Rudder Angle : ±20 deg. 
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Figure 4.10 Performance of Real System, Compensated, Complex Wave Input 

/; = 1.0, k | = 2.0. K, = 2.0 

Maximum Rudder Rate : ± 10 deg 4ec, Maximum Rudder Angle : ±20 deg. 
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Figure 4.1 1 Performance of Real System, Compensated, Complex Wave Input 

Z = 1 .0, k , = 6.0, K, = 6.0 

Maximum Rudder Rate : ± 10 deg/Sec, Maximum Rudder Angle : ±20 deg. 
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Figure 4.12 Performance of Real Svstem. Compensated, Complex Wave Input 

Z= 1.0, R, = 8.0, K, = S.O K h 

Maximum Rudder Rate : ± 10 degree, Maximum Rudder Angle : ±20 deg. 
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Figure 4.13 Maximum Rudder Rate vs Maximum Roll Anele 
Maximum Rudder Angle : 5 deg 
<*) = 0.5 rad/sec. 
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Figure 4.14 Maximum Rudder Rate vs Maximum Roll Angle, 
Maximum Rudder Angle = 10 dee, 'sec 
0) = 0.5 rad/sec. 
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Figure 4.15 Maximum Rudder Rate vs Maximum Roll Angle, 
Maximum Rudder Angle = 20 deg 
to = 0.5 rad sec. 
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Figure 4.16 Maximum Rudder Rate vs Maximum Roll Angle, 
Maximum Rudder Ancle = 30 dec 
(0 = 0.5 rad sec. 
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Figure 4.17 Maximum Rudder Rate vs Maximum Heading Angle, 
Maximum Rudder Angle = 5 deg 
0 ) = 0.5 rad; sec. 
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Figure 4.18 Maximum Rudder Rate vs Maximum Heading Angle, 
Maximum Rudder Angle = 10 deg 
co = 0.5 rad/ sec. 
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Figure 4.19 Maximum Rudder Rate vs Maximum Headine Ancle, 
Maximum Rudder Ancle = 20 deg 
(0 = 0.5 rad; sec. 



54 




Figure 4.20 Maximum Rudder Rate vs Maximum Heading Angle, 
Maximum Rudder Anale = 30 deg 
0> = 0.5 rad sec. 
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OPEN SYSTEM VS IDEAL SYSTEM 




Figure 4.21 Performance of Roll Reduction : Frequency vs Maximum Roll Angle 

Z= 1.0, K. = 4.0. K, = 4.0 y 6 

Open System vs IdeaFRudder. 
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Figure 4.22 Performance of Roll Reduction : Frequency vs Maximum Roll Ancle 
„ , ^ Z= 1.0, K. = 4.0, K, = 4.0 6 

Real Rudder : Maximum Rudder Rate = ± 5 deg/sec. 
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REfiL RUDDER 




Figure 4.23 Performance of Roll Reduction : Frequency vs Maximum Roll Angle 

Z = 1.0, K, = 4.0, K, = 4.0 
Real Rudder : Maximum Rudder Rate = ± 10 deg/sec. 
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ROLL-MAX (DEG) 



Figure 4.24 



REAL RUDDER 




Performance of Roll Reduction : Frequency vs Maximum Roll Angle 
Z = 1.0, K, = 4.0, K., = 4.0 
Real Rudder : Maximurfi Rudder Rate = ± 15 deg/ sec. 
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ROLL-MAX (DEG) 



RERL RUDDER 




Figure 4.25 



Performance of Roll Reduction : Frequency vs Maximum Roll Angle 
Z= 1.0, K. = 4.0, K, = 4.0 
Real Rudder : Maximum Rudder Rate = ± 20 deg/sec. 
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V. CONCLUSION AND RECOMMENDATION 



A. CONCLUSION 

From the results of the Rudder Roll Stabilization (RRS) computer simulations, 
the following conclusions have been reached : 

(1) Simple models can be derived to describe how wave forcing input and rudder 
angle affect the roll angle. 

(2) A properly designed roll controller will minimize the roll angle response and 
provide the desired heading angle as well. 

(3) This thesis has presented some practical aspects of a Rudder Roll Stabilization 
system. It has been shown that the steering machine may act as a severe 
limitation. When the rudder speed is sufficiently increased good roll reduction 
can be realized on existing ships. When it is possible, by means of a modified 
rudder construction, to increase the effect of the rudder on the roll without 
changing its effect on the heading an even better performance can be 
obtained. 

(4) It appears that both the rudder angle limit and the rudder speed limit have to 
be chosen carefully in order to realize a reasonable roll reduction. That is, the 
effectiveness of a RRS system is directly dependent upon available rudder 
moment. The available moment, in turn, depends on rudder rate. RRS 
system performance is rudder rate limited. Another parameter which is very 
important with respect to the roll reduction achieved by means of the rudder is 
the gain of roll and roll rate feedback (see Tables 3, 5 and 7). 

(5) The primary objective of the RRS has been achieved. That is, maximum roll 
reduction up to 55 percent has been measured with an ideal rudder (Z= 1.0, 
Kj = 4.0, K 2 = 4.0), and maximum roll has also been reduced to 47.5 percent 
with a reasonable real rudder (Z=1.0, Kj=4.0, K 2 = 4.0, 5 max =10 deg/sec, 
5 max = 20 deg), (see Tables 4 and 5). 



61 



B. RECOMMENDATION 

The results and conclusions of the RRS computer simulation lead to the 
following recommendation : 

(1) The results of this thesis demonstrate that roll stabilization using the rudder is 
feasible. Only one roll controller has been studied. It used roll angle and roll 
rate feedback. It is recommended that roll controllers including roll 
acceleration feedback should be studied. 

(2) In this thesis it has been assumed that cross coupling between the parameters 
of roll and yaw rate ( N<I> ) is zero. However, it is recommended that cross 
coupling between these parameters be considered in further studies on this 
topic. 

(3) The study does not consider wave input effects on the yawing moment. 
Further research should consider these effects. 

(4) Future research should investigate the use of advanced adaptive control and 
optimal estimation techniques to solve the problem of Rudder Roll 
Stabilization. 

(5) Future research should investigate the effects of various sea state conditions 
on the Rudder Roll Stabilization. 
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APPENDIX A 

CONCEPT BEHIND THE RUDDER-ROLL-STABILIZATION 



Roll stabilization is actually a secondary function of the rudder. Its primary 
purpose is to steer. However, as shown in Figure A.l, the rudder produces 
simultaneously a roll moment and a yaw moment which is needed to change ship 
course. 

Typical ship response periods to a roll moment are about 8 to 12 seconds, where 
as, typical response periods a yaw moment are about 30 to 35 seconds. 

The significant difference between these response periods allows the simultaneous 
superposition of yaw and roll control signals on the rudder without adversely affecting 
the response in either mode. 

The rudder moment is proportional to of the ship velocity and to the rudder 
angle. Using the rudder for roll stabilization involves opposing the wave induced roll 
moment by the rudder induced roll moment. 

Stabilization is achieved by adding to the steering control signal a roll control 
siiznal. 




Figure A.l Rudder Moments and Ship Response. 
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APPENDIX B 

THE ROLL MODEL WITHOUT COMPENSATOR 




U)\ 



Figure B.l Block Diagram Model of Ship's Dynamics with Autopilot 

(without compensator)." 

where 

Gj is the integration function (= 1 / S } . 

G-> is the ship's heading equation ( = 1/(St + 1)}. 

G 3 is a constant relating rudder angle to heading rate { = N§}. 



64 



G 4 is the rudder servo equation { = KRSER/(S + KRSER)}. 

G 5 is the autopilot equation { = KHP + KHD*S}. 

G 6 is the ship's roll equation { = <o n 2 /(s 2 + 2^co n s + co n 2 )}. 

G ? is a constant relating rudder angle to roll moment { = K§}. 

G g is a constant relating rudder angle to heading rate { = N<p}- 

G 9 is a constant relating heading rate to roll moment { = K f }. 

C, is a damping ratio of the roll dynamics. 
co n is the natural frequency of the roll dynamics. 

KRSER is a constant relating rudder servo gain. 

KHP is a constant relating heading gain (Autopilot). 

KHD is a constant relating heading rate gain (Autopilot). 



letting to., = 0.0, 

For convenience, products of the types GjG^G^... are written as Gj 2 3 ..., 



+ A - B) 



(B.l) 



A = - G 



1.2,4, 5, 7 



D 



(B.2) 




(B.3) 



C= - G 



1,2, 3, 4, 5 



D 



(B.4) 



D = C - E 



(B.5) 
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e=g 8 <» 



(B.6) 



Thus 



D= ■ G li2t3t4iS D - G 8 0 



-g 8 o> 



1 + G 



1,2, 3, 4, 5 



^ G 6 ^ C0 1 ‘ G 1 , 2 , 4 , 5 , 7 D ' G 2 , 9 D ) 



= GJ<0 + <£uA5isl2iML<l> ) 
61 1 1 + G 



1 . 2 , 3 , 4,5 



(B.7) 



(B.8) 



^_ = G 6< 1 + G 1 . 2 . 3 . 4 . 5 ) 

C0 1 1 + G 1 , 2 , 3 , 4,5 ‘ G l , 2 , 4 , 5 , 6 . 7, 8 “ G 2 , 6 , 8,9 

where G s is zero 



(B.9) 



<D 




(B.10) 



The physical interpretation of this is that in a beam sea, the waves do not cause 
yawing, and there is no cross coupling from roll to yaw, so the ship stays on course 
and the autopilot does not need to give rudder commands. Thus ship rolling is not 
affected by the rudder, because the rudder remains in the neutral position. 
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APPENDIX C 

THE ROLL MODEL WITH COMPENSATOR 




Figure C. 1 Block Diagram Model of Ship's Dynamics with Autopilot 

(With Compensator). 



where 

Gj is the integration function {= l/S}. 

G 2 is the ship's heading equation ( = l/(St + 1)}. 

G 3 is a constant relating rudder angle to heading rate { = Ng}. 
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G 4 is the rudder servo equation {= KRSER/(S + KRSER)}. 

G $ is the autopilot equation { = KHP+ KHD*S}. 

G 6 is the ship's roll equation ( = co n 2 /(s 2 + 2^co n s + co n 2 )}. 

G 7 is a constant relating rudder angle to roll moment (= Kg}. 

G g is a constant relating rudder angle to heading rate { = N ( j ) }. 

G 9 is a constant relating heading rate to roll moment ( = K f }. 

G c is a roll compensated equation (K 2 *S + Kj}. 

C, is a damping ratio of the roll dynamics. 
co n is the natural frequency of the roll dynamics. 

KRSER is a constant relating rudder servo gain. 

KHP is a constant relating heading gain (Autopilot). 

KHD is a constant relating heading rate gain (Autopilot). 



letting O)-, = 0.0, 

For convenience, products of the types GjG^^... are written as Gj 2 y.., 



<P = G 6 (co 1 + A - B) 



(C.l) 




(C.2) 




(C.3) 




(C.4) 
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D = C - E 



(C.5) 



e=g 8 o> 



(C.6) 



Thus 



D ' G 1,2,3,4,5 D ‘ G 3,4, G c <I> ' G 8 <1> 



- - ^ G 3.4. G c + G 8>* 
1 + G 1,2,3,4,5 



(C.7) 



<t> G 6 ( <0j - Gj 24 5 ? D - G 4 i 7 G c O - G 29 D) 



= G-[ 0). + t (G l,2.3,4,4,sA +G ^?,4,?,7,? + G ?,LiL22f. + - G 2,8,9> . G G ) <p , 

o l 1 i + rz 4,7, c J J 

1 U 1A 3,4,5 



0) 



G 6^ +G 1.2.3.4.5^ 



0) 1 1 + G 1,2,3,4,5 ' G l,2,4,5,6,7,8 ' G 2,6,8,9 + G c( G 4,6,7' G 2,3,4,6,9) 



(C.9) 



where G g is zero 



<P 



(0 



G 6^ 1 + G 1.2.3.4.S^ 

1 1+G 1,2,3,4,5 + G c< G 4,6,7 - G 2,3,4,6,9^ 



(C.10) 



(C.S) 
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APPENDIX D 

COMPUTER PROGRAM FOR ROLL, UNCOMPENSATED 



TITLE Model Without Compensator 

**************************************************************** 
*** *** 

*** The Roll Response of the Open System *** 

*** 'k'k’k 

******** ***** * ********************************************** **** 
* 

CONST ZETA=0 . 18 , 0MEGAN=0 . 59 , . . . 

KWH=0 . 0 , KWR=1 . 0 

* 

PARAM W=0 . 3 

PARAH A0=0 . 06 
* 

* 

INITIAL 

Y0=0.0 
YD0=0 . 0 

* 

* The Wave Equation 

* 

TH=W*TIME 
WAVE=A0*SIN (TH) 

* The Roll Equation 

ERR0R=(0MEGAN**2*WAVE*KWR)-SUM 

SUM= ( 2 *ZETA*0MEGAN*R0RATE ) + ( 0MEGAN**2*R0LL ) 

R0ACC=ERR0R 

RORATE=INTGRL (Y0 , ROACC) 
rorang=rorate*180 . /3 . 14 
R0LL=IMTGRL ( Y0 , RORATE ) 

ROLANG=ROLL*180./3.14 

CONTRL FINTIM=500 . 0 , DELT=0.1 
PRINT 1.0, rorang , ROLANG 

SAVE 0 . 1 , rorang, ROLANG 

* 

END 

PARAM W=0 . 4 

PARAM A0=0 . 10 

* 

END 

PARAM W=0 . 5 

PARAM A0=0 . 15 
* 

END 

PARAM W=0 . 6 

PARAM A0=0 . 10 
* 

END 

PARAM W=0 . 7 

PARAM A0=0 . 09 
* 

END 

PARAM W=0 . 8 

PARAM A0=0 . 07 

* 

END 

PARAM W=0 . 9 
PARAM A0=0 . 07 
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* 

END 

PARAM W=1.0 

PARAH A0=0 . 05 
* 

END 

PARAM W=l.l 

PARAM A0=0 . 04 
* 

* 

GRAPH (G,DE=TEK618,PO=0, NI=8) TIME(UN=SEC) ROLANG ( UN=DEG ) 

LABEL (G) ROLL ANGLE 

END 

STOP 
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